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We present measurements of the branching fractions for B° meson decays to rj' K° and uiK°, 
and of the branching fractions and CP- violation charge asymmetries for B + meson decays to r)-K + , 
r)K + , rj'Tr + , r]'K + , uitv + , and uK + . The data, collected with the BABAR detector at the Stanford 
Linear Accelerator Center, represent 383 million BB pairs produced in e + e~ annihilation. The 
measurements agree with previous results; we find no evidence for direct CP violation. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh, 12.39.St 



Charmless B decays are becoming increasingly useful 
to test the accuracy of theoretical estimation methods, 
such as those based on QCD factorization [l|, |2j, |3[ or fla- 
vor SU(3) symmetry [3,H, Q. In this paper we present 
measurements of branching fractions and, where appli- 
cable, charge asymmetries, for eight charmless B de- 
cays (and their charge-conjugates, implied throughout 
the paper): B+ -> t]tt+ , B+ -> r)K + , B+ -> 77'^+, 
B+ -► r)'K + , P° -► ri'K , B+ -► lott+ , B+ -► u>K+, 
and P° — > wiC . The results presented here represent 
improvement in precision overprevious measurements of 
these quantities by BABAR @,lEll, Belle [Iq|,EiLE1 and 
CLEO We previously reported a branching frac- 

tion limit for B° — > 77X [Tf, and CP asymmetries for 
B° -> 77'if and P° -> wif°ja[ii|. 

Charmless P decays with kaons are usually expected 
to be dominated by b — > s loop ("penguin") amplitudes, 
while 6 — > m tree amplitudes typically dominate for the 
decays with pions. However, the B — > rjK decays are 
especially interesting since they are suppressed relative 
to the abundant B — > rj'K decays due to destructive 
interference between two penguin amplitudes |T(> ■ The 
Cabibbo-Kobayashi-Maskawa (CKM) suppressed b — > u 
tree amplitudes may interfere significantly with 6 — > s 
penguin amplitudes of similar magnitudes, possibly lead- 
ing to large direct CP violation in B + — > 7771- + and 
P + — > ?y'7r + Jl7[i numerical estimates are available in 
a few cases [ITS; 0, El- We search for such direct 
CP violation by measuring the charge asymmetry A c h = 
(T- - r+)/(r- + T+) in the rates T± = T(B ± -> / ± ) 
for each charged final state / ± . 

Finally, phenomenological fits to the branching frac- 
tions and charge asymmetries of charmless B decays can 
be used to understand the relative importance of tree 
and penguin contributions and may provide sensitivity 
to the CKM angle 7 0, H, @, El , or to the effect of non- 
Standard-Model heavy particles in the loops [20| . 

The results presented here are based on data collected 
with the BABAR detector [5l| at the PEP-II e+e~ col- 
lider 22J located at the Stanford Linear Accelerator Cen- 
ter. An integrated luminosity of 347 fb , corresponding 
to 383 x I0 6 BB pairs, was recorded at the Y(4S) reso- 
nance (center-of-mass energy y/s = 10.58 GcV). 

Charged particles from the e + e~ interactions are de- 
tected, and their momenta measured, by a combination 
of five layers of double-sided silicon microstrip detectors 
and a 40-layer drift chamber, both operating in the 1.5 T 
magnetic field of a superconducting solenoid. Photons 



and electrons are identified with a CsI(Tl) electromag- 
netic calorimeter (EMC). Further charged particle iden- 
tification (PID) is provided by the average energy loss 
(dE/dx) in the tracking devices and by an internally re- 
flecting ring imaging Cherenkov detector (DIRC) cover- 
ing the central region. 

We establish the event selection criteria with the aid of 
a detailed Monte Carlo (MC) simulation of the B produc- 
tion and decay sequences, and of the detector response 
[23I ]. These criteria are designed to retain signal events 
with high efficiency. When applied to the data, they re- 
sult in a sample much larger than the expected signal, 
but with well characterized backgrounds. We extract the 
signal yields from this sample with a maximum likelihood 
(ML) fit. 

The P-daughter candidates are reconstructed through 
their decays ir° — > 77, K° — > K2 — * 7r + 7r~, u> — > 



(TEbr), 
o 



V 



7T + 7T 7T , T) — * 77 (?7 77 ), Tj — > TT + 7T~ 

?7 77 7r+7r- (rf n „„), and 77' -► p° 7 (77^), where p° 
The invariant mass of these particles' final states are re- 
quired to lie within about two standard deviations of the 
nominal mass 24] unless the mass is an observable in 



the ML fit, in which case we accept a wider range. For 
a Kg candidate we require a successful fit of the decay 
vertex with the flight direction constrained to the pion 
pair momentum direction, yielding a flight length greater 
than three times its uncertainty. Secondary charged pi- 
ons in 7/, r\ and to candidates are rejected if classified as 
protons, kaons, or electrons by their DIRC, dE/dx, and 
EMC PID signatures. For the primary charged track 
in B + decays we define the PID variables S„ and Sk as 
the number of standard deviations between the measured 
DIRC Cherenkov angle and that expected for pions and 
kaons, respectively. We include these observables in the 
ML fits to distinguish between primary 7r and K. For 
B + — > rf K + the backgrounds, including cross feed from 
the pion channel, are small. For this mode we perform a 
dedicated fit with less restrictive continuum background 
rejection (see below), and Sk < 2 to exclude pions (and 
lighter particles). 

We reconstruct the P-meson candidate by combining 
the four-momenta of a pair of daughter mesons, with a 
vertex constraint if the ultimate final state includes at 
least two charged particles. Since the natural widths of 
the 77, r;', and ir° are much smaller than the resolution, 
we also constrain their masses to nominal values [24| in 
the fit of the B candidate. From the kinematics of T(4S) 
decay we determine the energy-substituted mass tties = 
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-\J\s — p B and energy difference AE — E B — \\/~s, where 
(Eb,Pb) is the B-meson 4-momentum vector, and aii 
values are expressed in the T(4S') frame. The resolution 
in toes is 3.0 MeV and in AE is 24-50 MeV, depending 
on the decay mode. We require 5.25 < toes < 5.29 GeV 
and \AE\ < 0.2 GeV. 

Backgrounds arise primarily from random combina- 
tions of particles in continuum e + e~ — > qq events [q = 
u,d,s,c). We reduce these with requirements on the an- 
gle 9t between the thrust axis of the B candidate in the 
T(4S) frame and that of the rest of the charged tracks 
and neutral calorimeter clusters in the event. The dis- 
tribution is sharply peaked near | cos 6* T | = 1 for qq jet 
pairs, and nearly uniform for -B-meson decays. We re- 
quire | cos6» T | < 0.90 (< 0.65 for r)' p ^Tr+, < 0.80 for lott+ 
and u>K + ), which optimizes the expected signal yield rel- 
ative to its background-dominated statistical error. In 
the ML fit we discriminate further against qq background 
with a Fisher discriminant T that combines several vari- 
ables which characterize the energy flow in the event [25| . 
It provides about one standard deviation of separation 
between B decay events and qq background (see Fig. [lb). 

We also impose restrictions on resonance decay an- 
gles to exclude the most asymmetric decays where soft- 
particle backgrounds accumulate and the acceptance 
changes rapidly. We define the decay angle 9 r Acc for a 
meson r that decays to two particles as the angle be- 
tween the momenta of a daughter particle and the me- 
son's parent, measured in the meson's rest frame. We de- 
fine W = cos8 r dcc and require \H P "\ < 0.9 for B -> rf^K 

and \n p °\ < 0.7 for B+ v' P1 ^ + ■ For the three-body 
lu — > 3-7T mode the direction for the decay is the normal 
to the decay plane, and we include as an observable 
in the ML fit. 

The average number of candidates found per selected 
event is in the range 1.05 to 1.13, depending on the fi- 
nal state. We choose the candidate with the daughter 
resonance mass closest to the nominal value. From the 
simulation we find that this algorithm selects the correct- 
combination candidate in about two thirds of the events 
containing multiple candidates, and that it induces neg- 
ligible bias in the ML fits. 

We obtain yields for each channel from an extended 
maximum likelihood fit with the input observables AE, 
™ES: F , m r (the invariant mass of the rj, r( , or lu candi- 
date), and, for charged decays other than B + — > r]'K + , 
the PID variables «SV and Sk- The selected data sample 
sizes are given in the second column of Table QJ Besides 
the signal events they contain qq (dominant) and BB 
with b — *■ c combinatorial background, and a fraction of 
background from other charmless BB modes, which we 
estimate from the simulation to be less than 2% of the to- 
tal fit sample. The latter events have ultimate final states 
different from the signal, but with similar kinematics so 
that broad peaks near those of the signal appear in some 



observables, requiring a separate component in the prob- 
ability density function (PDF). The yield of this compo- 
nent is free in the fit for all cases except B° — ► loK®, 
where the fit stability requires fixing the yield to the ex- 
pectation from MC. The likelihood function is 

N 

£ = exp(-E^)lIE^ x 

j,k i j,k 

V^m^VjinV^AEl) [V^SDV^VM)] , 

where N is the number of events in the sample, and for 
each component j (signal, combinatorial background, or 
charmless BB background) and flavor k (primary K + or 
7T+), Yjk is the yield of events and Vj {x].) the PDF for ob- 
servable Xk in event i. Some factors in [ ] are omitted for 
some modes. The flavor-dependent factors T'j{AE l k ) and 
Vj(S l k ) take common functional forms for pion or kaon, 
e.g., Fj(AEi) or Fj (AE^ = AE^ + SAE{p 1 )), where p 
is the primary track momentum; SI is treated similarly. 
For the modes B — » rj'^K we found no need for the BB 
background component. The factored form of the PDF 
indicated in Eq. [1] is a good approximation, particularly 
for the qq component, since correlations among observ- 
ables measured in the data are typically a few percent or 
less. Distortions of the fit results caused by our approxi- 
mations are measured in simulation and included in the 
bias corrections and systematic errors discussed below. 

We determine the PDFs for the signal and BB back- 
ground components from fits to MC samples. We cali- 
brate the resolutions in AE and toes with large data con- 
trol samples of B decays to charmed final states of simi- 
lar topology (e.g. B — > D{K-k-k)-k). We develop PDFs for 
the combinatorial background with fits to the data from 
which the signal region (5.27 GeV < toes < 5.29 GeV 
and \AE\ < 0.1 GeV) has been excluded. 

We use the following functional forms for the PDFs: 
sum of two Gaussians for V a i g (m-Es), BB (AE), and 
the sharper structures in V BB (rriEs) and Vj{rn r )\ linear 
or quadratic dependences for combinatorial components 
of V BB q q(m r ) and for VqqiAE); and a Gaussian func- 
tion with separate low- and high-side width parameters 
for Vj(F). The qq background in toes is described by 
the threshold function x\J\ — x 2 exp [— £(1 — £ 2 )], with 
x = 2toes/a/s and parameter £. These functions are 
discussed in more detail in [25], and some of them are 
illustrated in Fig. [T] 

We allow the parameters most important for the de- 
termination of the background PDFs to vary in the fit, 
along with the yields for all components, and for charged 
modes the signal and qq background charge asymmetries. 
Specifically, the free background parameters are most or 
all of the following, depending on the decay mode: £ 
for toes, linear and quadratic coefficients for AE, area 
and slope of the combinatorial component for m r , and 
the mean, width, and width difference parameters for T . 
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Results for the signal yields are presented in the third 
column of Table [T] for each sample. 

We validate the fitting procedure by applying it to en- 
sembles of simulated qq experiments drawn from the PDF 
into which we have embedded the expected number of sig- 
nal and BB background events randomly extracted from 
the fully simulated MC samples. Biases obtained by this 
procedure with inputs that reproduce the yields found in 
the data are reported in the fourth column of Table |TJ 

In Fig. Q] we show, as a representative of the fits, the 
projections of the PDF and data for the B + — ► uK + fit, 
and in Fig. [2] projections onto rriES for each of the eight 
decays, with submodes combined. The data plotted are 
subsamples enriched in signal with a threshold require- 
ment on the ratio of signal to total likelihood (computed 
without the plotted variable) that retains 35%-80% of 
the signal, depending on the mode. 

We determine the reconstruction efficiencies as the ra- 
tio of reconstructed and accepted events in simulation to 
the number generated. We compute the branching frac- 
tion for each channel by subtracting the fit bias from the 
measured yield, and dividing the result by the efficiency 
(including secondary branching fractions) and the num- 
ber of produced BB pairs [25|. We assume equal decay 
rates of the T(AS) to B + B~ and B°B° . Table [Q gives the 
numbers pertinent to these computations. The statistical 
error on the signal yield or branching fraction is taken as 
the change in the central value when the quantity — 2 InC 
increases by one unit from its minimum value. 

We combine results where we have multi- 
ple decay channels by adding the functions 
-21n{[£(B)/£(B )]®G(cr')}, where B is the cen- 
tral value from the fit for each decay channel, and <S)G 
denotes convolution with a Gaussian function to include 
the systematic error a' discussed below. We give the 
resulting final branching fractions for each mode in 
Table H 

Systematic uncertainties on the branching fractions 
arise from the PDFs, BB backgrounds, fit bias, and effi- 
ciency. PDF uncertainties not already accounted for by 
free parameters in the fit are estimated from the consis- 
tency of fits to MC and data in control modes. Varying 
the signal-PDF parameters within these errors, we es- 
timate yield uncertainties of 0.4%-2.2%, depending on 
the mode. For the BB backgrounds we vary the in- 
put branching fractions within their uncertainties for the 
modes that contribute most to the selected sample. The 
resulting changes in the signal yield are taken in quadra- 
ture and scaled to the total of all modes to determine the 
systematic uncertainty. For the r/'^K modes, where no 
BB component is used, we use 10% of the expected BB 
background in the sample as this is the typical correla- 
tion with the signal yield. For ujK® where the BB yield 
is fixed, we take as a systematic uncertainty the average 
change in the signal yield when the BB yield is varied be- 
tween zero and twice the nominal value. The uncertainty 
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FIG. 1: Plots of signal-enhanced subsets of the data distribu- 
tion for B + — > wi^projected on each of the fit variables: (a) 
«ies, (b) A_E, (c) T, (d) "H° , (e) u> mass, and (f) St- Points 
with errors represent the data, solid curves the full fit func- 
tions, dashed curves the sum of the background functions, and 
dot-dashed curves the signal from B + — » utv + . The variable 
AE is computed with the pion mass. 



of the bias (Table Q} is a quadrature sum of its com- 
ponents: the statistical uncertainty from the simulated 
experiments, and half of the corrections attributable to 
correlations omitted from the signal and BB background 
models, and to PID of the primary charged track. The 
primary-track PID correction is significant only for mis- 



identified kaons from B + — ► rj'K + in the B + 
channels. 

Uncertainties in our knowledge of the efficiency, found 
from auxiliary studies, include 0.5% x 7V t and 1.5% x 7V 7 , 
where N t and 7V 7 are the number of tracks and photons, 
respectively, in the B candidate. The uncertainty in the 
total number of BB pairs in the data sample is 1.1%. 
Published data 24| provide the uncertainties in the B- 
daughter product branching fractions (0.7-3.2%). The 
uncertainties in the efficiency from the event selection 
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TABLE I: Number of events N in the sample, fitted signal yield Is, and measured bias (to be subtracted from Ys) in 
events (ev.), detection efficiency e, daughter branching fraction product (n^Oi an d measured branching fraction B and charge 
asymmetry Ach with statistical error for each decay chain, and for the combined measurements the branching fraction and 
charge asymmetry with statistical and systematic error. The number of produced BB pairs is given in the text. 



Mode 


AT fair \ 


is (ev.j 


Bias (ev. 


e (%) na (%) 


k? /-in— 6\ 
(10 j 


A 

■Ach 














5.0 ± 0.5 ± 0.3 


-0.08 ± 0.10 ± 0.01 




44883 


oco+30 
ZOo 29 


6±3 


04;.! 




4.9 ±0.6 


-0.05 ±0.12 




22333 




6±3 


ZO.O 


oo a 


5.5 ± 1.0 


-0.13 ±0.18 


T)K + 












3.7 ± 0.4 ± 0.1 


-0.22 ± 0.11 ± 0.01 


T} 11 K + 


44883 


197+11 


6±3 




QO A 


3.9 ±0.5 


-0.25 ±0.13 




22333 


7i±iS 


4±2 


9Q 9 




3.3 ±0.8 


-0.15 ±0.23 


77' 71+ 












3.9 ± 0.7 ± 0.3 


0.21 ± 0.17 ± 0.01 




16879 




14±3 


27.2 


17.5 


4.0 ± 0.9 


0.14 ± 0.20 


77p 7 7T + 


35523 


97tg 


23±7 


18.4 


29.4 


3.6 ± 1.1 


0.35 ±0.30 


■n'K+ 












70.0 ± 1.5 ± 2.8 


0.010 ± 0.022 ± 0.006 


ri K+ 

' 1 7] TV TV 1 v 


3170 


1060±35 


Oil 


23.2 


17.5 


68.2 ±2.3 


-0.005 ± 0.033 




79501 


2405±69 


31±16 


29.2 


29.4 


72.2 ±2.1 


0.022 ± 0.028 


r,'K° 












66.6 ± 2.6 ± 2.8 


(see [15]) 




1100 


329±20 


3±1 


23.2 


6.1 


60.7 ±3.7 




T)'K° 


19927 


831±38 


35±17 


28.0 


10.2 


72.8 ±3.5 




tJ7T+ 


76735 


516±38 


44±22 


20.5 


89.1 


6.7 ± 0.5 ± 0.4 


-0.02 ± 0.08 ± 0.01 


u:K+ 


76735 


457±32 


29±15 


20.0 


89.1 


6.3 ± 0.5 ± 0.3 


-0.01 ± 0.07 ± 0.01 




15914 


146±18 


10±5 


21.2 


30.8 


5.4 ± 0.8 ± 0.3 


(see [9J) 




k25 5.26 5.27 5.28 5.29 
m ES (GeV) 



£25 5.26 5.27 5.28 5.29 
m ES (GeV) 



the interactions with material in the detector. We study 
these by comparing this effect in MC for the signal, qq 
background in the data, and control samples mentioned 
previously. We apply corrections, and assign systematic 
errors, to A c h equal to —0.010 ± 0.005 for modes with a 
primary kaon and 0.000 ± 0.005 for those with a primary 
pion. We apply an additional correction with uncertainty 
for dilution of the A c h measurement associated with the 
yield bias, which is significant only for B + — > rj' n ir + . 
This is obtained from the same MC studies that are used 
to estimate the yield bias. 



After combining the measurements we obtain for the 
branching fractions: 



FIG. 2: Plots of signal-enhanced subsets of the data distribu- 
tions projected onto ?ties for the decays: (a) B + — > 7771+, 
(b) B+ -» rjK+, (c) B+ -> i 1 '-k+, (d) B+ -> 77' K + , (e) 
5° -> rj'K", (f) B+ -> U7V+, (g) B + -> loK+, and (h) 
B° — » luK°. The solid line represents the result of the fit, 
and the dot-dashed line the background contribution. The 
dashed line gives the sum of background and the 773^ (a, b) 
or ri' V7rTr (c-e) component of the signal. The dotted line shows 
the K or n cross-feed component, where applicable. 



are below 0.5%. 

For the measurements of A c h , biases arise in principle 
from charge-dependent effects in the track reconstruction 
or particle identification, or from imperfect modeling of 



B(B+ 


-> 7/7T+) = 


(5.0 ±0.5 ±0.3) x 10~ 6 


B{B+ - 


■+ rjK+) = 


(3.7 ±0.4 ±0.1) x 10~ 6 


B{B+ - 


->77'7T + ) = 


(3.9 ±0.7 ±0.3) x 10~ 6 


B{B+ - 


>r}'K + ) = 


(70.0 ±1.5 ±2.8) x 10~ 6 


B(B° - 


♦ V 'K°) = 


(66.6 ±2.6 ±2.8) x 10~ 6 


B{B+ - 


^W7T+) = 


(6.7 ±0.5 ±0.4) x 10~ 6 


B(B+ - 


>LUK+) = 


(6.3 ±0.5 ±0.3) x 10~ 6 


B(B°- 


-+u;K ) = 


(5.4 ±0.8 ±0.3) x 10" 6 . 
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For the charge asymmetries we find 



A ch {B + 


— > rjn + ) 


= -0.08 ±0.10 ±0.01 


A ch {B+ - 


-4 


= -0.22 ±0.11 ±0.01 


A ch (B+ - 


-> r/'7r+) 


= 0.21 ±0.17 ±0.01 


A ch (B+ - 


* T)'K + ) 


= 0.010 ±0.022 ±0.006 


A ch {B+ 


—> UJTT + ) 


= -0.02 ±0.08 ±0.01 


A ch (B+ - 


->luK+) 


= -0.01 ±0.07 ±0.01. 



The first error quoted is statistical and the second sys- 
tematic. These results are generally consistent with pub- 
lished measurements Rfl S3, Ell Ell [3 Ell] an d supersede 
our previous ones EL la HI; for B(B + — > r]K + ) we find a 
value about twice that of [lOj. The theoretical estimates 
are in agreement with the data (though the data have 
been used in some predictions), but with greater uncer- 
tainty [3, 0,0]. Approaches that fit all available data with 
a moderate number of model parameters have proved 
fruitful [3, S We find no clear evidence for direct 
CP-violation charge asymmetries in these decays. The 
world average of the measurements of A c n for B + — > r/7r + 
(B + — > r]K + ) are both negative and 2.3 (3.0) standard 
deviations from zero, while the predictions of 0] are pos- 
itive, though with large errors. 
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